Quasilinear diffusion for the chaotic motion of a particle in a set of longitudinal waves* 
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The rigorous analytical calculation of the diffusion coeffi- 
cient is performed for the chaotic motion of a particle in a set 
of longitudinal waves with random phases and large ampli- 
tudes (~ yl). A first step proves the existence of a quasilinear 
diffusion on a time scale ~ A~'^^^ In A. A second step uses 
this property to extend the result to asymptotic times by in- 
troducing the conditional probability distribution of position 
and velocity of an orbit at a given time when they are known 
at a previous time. 
PACS numbers: 

52.35.-g (Waves, oscillations, and instabilities in plasmas and 
intense beams) 

05.45.-a (Nonlinear dynamics and nonlinear dynamical sys- 
tems) 

05.60.-k (Transport processes) 

52.20.-j (Elementary processes in plasmas) 



I. INTRODUCTION 

Many chaotic Hamiltonian systems encountered in 
physics display a chaotic diffusion and in many cases the 
corresponding diffusion coefficient is given by a so-called 
quasilinear estimate [0-1). The proof that this estimate 
is correct exists for the standard map with large con- 
trol parameter , but is lacking for other systems with 
a spatially smooth force. We provide this proof for the 
one-dimensional chaotic motion of a particle in a general 
set of waves. 

This result strengthens the link between the micro- 
scopic deterministic (chaotic) dynamics and the macro- 
scopic stochastic motion. Its extension to the self- 
consistent many-body problem is a central problem to 
non-equilibrium statistical physics. 

This paper is organized as follows. We first introduce 
our model dynamics and stress the core of our argu- 
ment. Then we recall the traditional argument deriv- 
ing the quasilinear diffusion over a time short with re- 
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spect to a characteristic time Tsprcad ~ A"'^/^ {A being a 
typical value of the wave amplitudes) and introduce the 
explicit form of the quasilinear diffusion coefficient. We 
rederive this result within our new approach and take ad- 
vantage of a better understanding of the particle motion 
to extend the validity of quasilinear diffusion to a time 
scale ^ A~'^/'^\!iA, which is longer than the traditional 
scale Tsproad for A large. Finally, we introduce the con- 
ditional probability distribution of position and velocity 
of a chaotic orbit at a given time when they are known 
at a previous time and, thanks to the non-confinement 
of the velocity of the chaotic orbit, we further extend the 
quasilinear estimate to asymptotic time scales. 

II. DYNAMICAL MODEL AND ASSUMPTIONS 

We consider the dynamics of a particle in a set of longi- 
tudinal waves (e.g. Langmuir waves) with random phases 
and large amplitude, as defined by the Hamiltonian 



H{p, q,t) 
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M 



^ ArnCOs{krnq - UJmt + ipm), (1) 



where the ifm^s are random variables, and the 
uJmYs are prescribed triplets of positive param- 
eters. Such a dynamical system has already been studied 
in the literature, and for large A^'s the diffusion coef- 
ficient has been found numerically to take on the quasi- 
linear value |,|-§ defined below [0. The average over 
M '3> A^^^lnA ^ 1 random phases is central to our 
proof, in agreement with the occurrence of uncontrolled 
phases in many experiments and with the fact that the 
transport in (|l|) is much less diffusion-like if one aver- 
ages only over initial conditions (poi'Zo) The large 
A limit (dynamically speaking, the limit of strong res- 
onance overlap parameter) corresponds to the limit of 
continuous spectrum often encountered in physics. 

In agreement with most of the literature on quasilin- 
ear transport, the analysis is performed here in terms 
of quadratic means, and not in terms of the probabil- 
ity distribution functions, but we indicate at the end of 
this paper how our technique could be used to prove the 
gaussianity of such functions. 

The equations of motion are 
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We first consider the time to be short enough for the orbit 
to stay close to the unperturbed orbit g^"-* (t) = go + Pot, 
and let Aq{t) = q{t) - q^°\t), Ap{t) = p(t) - pq. Wc 
compute their statistical properties when averaging over 
ah (/3„'s. 

For completeness, we first evaluate Ap{t) by the tra- 
ditional technique using first order perturbation in 
the amplitudes: 



ra—1 

[C0s{kmqo + (Pm) " COs(f2mi + kmqo + (frn)], (4) 

with f2„j — kmPo — oJm] if — for some to, 

the corresponding term in the sum is the limit as 
0. At this order, {Ap{t)) = and {Ap'^{t)) = 

Yj^=li^rnhnl^mf[l " C0s(f2„^t)]. 

Let Vm = ^^m/km- Wc assume that Ail^ — r2^_}_i — 
and Avm — I'm+i — Vm have a sign independent of to, 
which is natural for Langmuir waves and for the dynamics 
of Ref. a. Let 
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Dm may fluctuate with to, but we assume (for simplicity 
only) that for 



^ > 0, T.-=-L D 



I ^'^m-\-j I 



-L+1-Vm-L \ IS a 



constant -Dql, called the quasilinear diffusion coefficient. 
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Let AVLlm = max|rim+L+i 

and Tc (fimax - ^^min)^^ ; Tdiscr and Tc are respectively 
the discretization time and the correlation time of the 
wave spectrum as seen by the particle. 



III. NON-CHAOTIC INITIAL QUASILINEAR 
TRANSPORT 

Assuming Tc <C t <C Tjiscr, one obtains (Ap(<)^) = 
{2Dqi^/tt) - cos{Vtt)]<m = 2DQLi, where the 

discrete sum has been turned into an integral. As a 
result, the diffusion coefficient takes on the quasilinear 
value 1?QL- A similar calculation for q yields {Aq{t)) = 
and {Aq{tf) = 2DQLi^/3. For t < Tc, Ap grows 
linearly with time, and (Ap^) grows quadratically, as 
all modes act with a constant force on the orbit. For 
Tc <C t <C Tdiscr, the range of to contributing to the dif- 
fusion (modes acting with a nearly constant force) nar- 
rows like 1/t. The range of t is further restricted by 
the condition for the orbit to remain close to the unper- 
turbed one. This is traditionally obtained by requiring 
(kl^^^Aq^it)) < 47r2, namely t < T^proad with 



Tsprcad 



1/3 



(6) 



where we introduce the resonance broadening frequency 
7Dri = {k'^Dqi,)^/^ and take 7d = max„7D„. 

In our approach, we evaluate Ap{t) as in Ref. by 
integrating formally the equation of motion for p. This 
yields {Ap{t)) = over the range < i <C tql defined 
below, and {Ap'^{t)) = Ao + A+ + A_, with 
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where $,n(i) ~ k„iAq{t) -I- r2,„t + (^„i, with ri± = ±1 
and rjQ = —1, and under condition toi 7^ TO2 for 
J = — , and condition toi = m2 for j = 0. Let 
t_ = <i - i2 and t+ = (<i + i2)/2. For <_ < Tgproad, 
(exp[ifc„(Ag(t+ +t_/2) - A9(i+ -t_/2))]) may be con- 
sidered as equal to 1. Therefore the support in t_ 
of the integrand in Ag is of the order of Tc. We as- 
sume Tc <C Tsprcad- Hence the integration domain in 
t_ may be restricted to \t- \ < vTc where v is &. few 
units. In the limit where vTc <^ t <ti Tdiscr, we obtain 
Ao = E™=i/o(2^™A)/o'''"(cos[ant-])Ar!„,di_di+ = 
2-Dql I]m=i('^^"i)~"^(siii[^^mi^Tc])Af2„it = 2DQi,t, with 
the discrete sum over m approximated by an integral. 

For t Tspread WC approximate q(t) by its unper- 
turbed value q^^\t). As this orbit does not depend 
on the phases, the averaged cosines in (|^) are zero for 
j — ±, and so are the A±'s. Then our second ap- 
proach shows again that the diffusion coefficient takes 
on the quasilinear value. {Aq'^{t)) too may be com- 
puted by integrating the equation of motion [ p^ . This 
involves calculating {Ap(ti)Ap(t2)) , in the same way as 
{Ap^{t)), and one recovers the traditional estimate for 
{Aq^{t)). This provides a way for introducing the condi- 
tion t <C Tsprcad without rcsorting to the traditional per- 
turbative approach, and shows that the usual quasilinear 
diffusion coefficient may be recovered independently by 
our second approach. 



IV. CHAOTIC TRAJECTORY SPREADING 

In fact our second approach is much more powerful. 
As was pointed out in Ref. ||], A± vanishes provided 
that the dependence of Ag over any N^p — 2 phases with 
all other phases fixed is weak, a condition far less strin- 
gent than the previous condition = M which led to 
i ^ Tsprcad- Reference Q estimated the upper bound in 
time of the initial quasilinear diffusion through numerical 
calculations for moderate values of the waves amplitude. 
Here we derive such a bound analytically for large enough 
amplitudes. 
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We measure these amplitudes by the parameter En = 
[2£>QLfc„| AD„|/7r]^/^ which corresponds to the typical 
electric field of a wave. A related dimensionless quantity 
characterizes our scaling, namely the Chirikov resonance 
overlap parameter 

sK) = 2K/2+4y+i]/|A«„| (8) 

or equivalently the ratio 

SK) = fc„|Av„|/7D„~5s-4/3 

of the frequency mismatch between neighbouring waves 
(in the frame of either wave) to their resonance broad- 
ening frequency. As these quantities depend on n, they 
characterize the dynamics locally. In the following, we 
are interested in the dense spectrum, or strong overlap, 
or large amplitude limit. To ensure a genuine scaling, we 
consider families of dynamics where — Ean and 
the reference amplitudes a„ are constant while E' — > cx3, 
or B{vn) = Bbn and the coefficients 6„ are constant while 

Apart from the small dimensionless parameter B, we 
also introduce the Kubo number JCc = Tc/rgpicad- The 
wide velocity spectrum of the waves ensures that /Cc ^ 1- 

The limit of interest is the joint limit /Cc — > and 
S — > (or /Cc ^ and s oo). 

A. Spreading due to a single random phase 

In order to avoid too heavy formulas, we give the ex- 
plicit derivation for the spreading due to one phase, and 
extend the result to two phases afterwards. To estimate 
this spreading we study how the orbit which is at {qo,Po) 
at i = is modified when phase ipn changes from to a 
finite value. Let {q^{t) ^ p^{t)) be the orbit for (/3„ = 0, let 
5qn{t) = q{t) ~ q,/it) and Spn{t) = Sqn{t) = p{t) - 
We assume t to be small enough so that kinax\Sqn{t)\ ^ tt. 
As Sqn{t) is small, we may linearize the motion 

Spnit) ~ F{t)5qn{t) + ^„/c„(sin*„(i) - sm^no{t)) 

(10) 

where F{t) = YlZ=i Am cos ^„^{t), with = 
kmq^{t) - LOmt + ipm and = knqr^it) - ujnt. Then 
( p^ ) and initial conditions ((5q„(0), (5p„(0)) = (0, 0) imply 

Sqnit) = f\t - t")F{t")5qn{t")dt" + %,o(i) (11) 
Jo 

where 

5g„o(i) = A„fc„/o/j'(sin*„(t")-sin*„o(i"))di"dt'. In 
the short-time limit, the dominant term in expression 
( |rT| ) for 5qn will be 5qnQ, but over longer times the first 
term may self-amplify and overtake the second one. 



We only estimate ((5q'„(t)^), but {Sqn{t)) can be com- 
puted by the same technique and turns out to be negli- 
gible over the time interval of interest. In a first stage, 
consider the contribution of Sqno to the variance, Co{t) ~ 
{Sq,,om = (A^A2/2)/„7„*'^/J/„*^(cos(vI/„(t'/) - 
'i'n{t2)))dt2dt'2dt'{dt'i. To estimate this expression, note 
that ^i!„{t'{)~^n{t'i)-n.n{t'{-t'^) = kn{qM)-Q¥it2))- 
knPoit'i ~ ^2)1 and, for the range of time of interest, 
qi/{t") ^ q'-^^t") is essentially the sum of M — 1 terms 
in which a random phase fm i'm- 7^ n) is added to a 
term which has a weak dependence on ipm- Therefore, 
this sum is almost gaussian, and for M ^ 1 we may ap- 
proximate q^(t") by a brownian motion. Furthermore, as 
M ^ 1, we approximate qyi{t") by q{t") in the averages. 
Using the distribution of Ag(t2) — Aq{ti), we find |1^] 
the estimate 

Coit) < CoM(t) = 0.28fc„|Ai;„|7g„t3 = 0.286 [jDntf. 

(12) 

For the second stage, we take into account the first 
term in the right hand side of ([lO|). As Sqn is small, 
we may treat F(t) as a gaussian process with moments 
(Fit)) = and {F{h)Fit2)) = 2jl,J{h - ^2) where 
S{t) is the Dirac distribution. Indeed q^{t) has a weak 
dependence on any phase (pm, which makes {F{ti)F{t2)) 
a Bragg-like function with the small width Tc in ti — 
t2. Higher moments of F are assumed to factorize, i.e. 
F is treated as a white noise, which is consistent with 
approximating q^(t) by a brownian motion. 

We estimate the spreading of Sqn{t) by computing 

C{t) EE {Sqnitf) 

^ rt rt, rt f\F[t'i)FH^i)){5qn{t'l)5qn{t'^)) 
Jo Jo Jo Jo 

dia'diadti'di'i + Co(t) 

= {Ey2) / / / C(i")dt"dt^dt'i 
Jo Jo Jo 

+Co{t). (13) 

It follows from (|o|) and our assumptions on F that 
C{t) ^ Coit)+LC{t) with 

Lfit) = (£;V2)/„7„77"(*-*^'/(i")dt"dt'2dfi. 

As {l-Ly^ preserves positivity ||l|l, C = {l-L)'^^Co < 
(1 — L)^^CoM = Cm- Applying the Laplace transform to 
both sides of equation Cm — Com + LCm, we compute 
Cm and find 

C(i) < CM(t) = O.UBk-^{e'' - 1 + 2g{t')) (14) 

with t' = A^/^'jBnt and g{t') = e-*'/^ cos(t'V3/2) - 1. 
This estimate for the variance of 6qn{t) starts from zero 
at t = and diverges exponentially for i — > 00. Its expo- 
nentiation time scale TLiap ~ ^ Tsprcad is the recip- 
rocal of the Liapunov characteristic instability rate (this 
is reminiscent of Ref. 0]). However, as the coefficient in 
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front of the exponential goes to zero as E —>■ oo, the time 
needed by our upper estimate on k'^C(t) to reach unity 
is of the order of 



^QL = 7d 



IhiKI 



(15) 



Though this time goes to zero as i? ^ oo, it is 0{\nB~^) 
times larger than the time Tspioad over which the initial 
quasilinear approximation is traditionally justified. 



B. Spreading due to two random phases 

The result of this discussion is that "(/(i) depends little 
on any given phase over a time tql". For Af ^ 1, the 
argument is easily strengthened into "g(i) depends little 
on any two given phases over a time tql". To this end 
{lvh.^,iiA^),Piit,,/i2{i)) and {5q{t), 5p{t)) are defined start- 
ing from ipmi = </'m2 = 0, and a third term similar to the 
second one adds in the right hand side of (|l^). The first 
stage of our iteration procedure now estimates the con- 
tribution of both phases and ^Pm2 by a term again 
of the order of B~'^^^t^ , while the second stage does not 
change. 

As a result, for t <C tql, the non-quasilinear terms 
A± are negligible since q has a small dependence on any 
given pair of phases in this time range. Furthermore these 
terms may be estimated by expliciting in the argument of 
the cosine of (|7|) the main dependence over and ipm2 
through estimates and <5<i>m2 of the type kmSqnO for 

both phases, and by expanding to second order in these 
(5<I>'s. Such estimates hold for t <C /3tql with < (3 < 1 
for E large enough, and yield A_|_ ~ E'^t^P and A_ ~ 
EH^ which are negligible with respect to Aq in the time 
interval of interest, and do not grow with M although 
there are 2Af ^ — M "off-diagonal" terms. 



V. QUASILINEAR TRANSPORT OVER LARGE 
TIMES 

Finally, we show that the quasilinear estimate holds 
for asymptotic times. Let Pmin = niin(t;m) and Pmax = 
max(wm). We assume that in the velocity domain 
[Pmin , Pmax] the dynamics is chaotic enough for a typi- 
cal orbit to be unconfined in p within this domain, but 
that the time of interest is also smaller than the time for 
the orbit to reach the boundaries of the chaotic domain. 
Therefore we set the condition min[(po — Pmin)^, [po — 

Pmax)^] > i^QLTQL ^ K'^lhn ln(-B"l) tO COmputC UOW 

the diffusion coefhcient due to the chaotic motion when 
M and E are large. We define 5q{T\p, q, t) = q{t + T) — q — 
pr, where q{t') is the position at time t' of an orbit which 
is at {p,q) at time t : 6q{T\p, q,t) tells the departure of 
this orbit from the free motion during the time interval 



Integrating formally the equation of motion for p yields 
(Ap^(t)) = - E E ^ ^"V"^" /* /*(cos$)dt'di" 

m,„=le=±l ^ -^0 Jo 

(16) 

where* = ik^+ek„)qit") + k„M[t' -t"\p{t"), q{t"),t"] + 
kmp{t"){t' — t")—uJmt' — ei^nt" + 'fm + i'Pn- We introducc 
the probability distribution P{Sp, i|po) of Sp = p{t) — po 
for an orbit started at p = po at to = 0; it is independent 
of go- 

{cos[k,n6q{t' - t"\p{t"),q{t"),t")]) is independent of 
q(t"), and its contribution for diagonal (m = n, e = —1) 
terms to (nm is 



B 



M 



hm y 



(^m ^m) 



2 /•* /•* 



P{Sp,t"\po) 



4t JO JO 
(cos[kmSq{t' - <"|po, q{t"), i")])* 
cos[fc™(po + 5p){t' - t") - uj„,{t' - t")] 
dSpAt'dt" 



M 



= lim > 



exp [iQmT] (exp [ifc„ Sq{T\po,q{t"),t")]), drdt" , (17) 

where the starred average means the average done with 
the constraint p{t") = Pq + Sp, and where the Fourier 
transform 



P{a,t"\po) 



P{Sp,t"\po)exp{ia5p)d6p (18) 



was used. As Sq is computed with the knowledge of p at 
time t" which sets only one condition on a set of many 
phases, an average with the constraint p{t") ~ po + Sp 
may be computed by using the initial quasilinear es- 
timate at time \t' — t"\ < tql. Hence the function 
{exp[ikfnSq{t' — t"\p, q{t"),t")])^ is correctly computed by 
the previous quasilinear estimate over its whole support 
in t' — t" as Tql 3> Tsprcad- This estimate is indepen- 
dent of p, and we could set p — po in the average co- 
sine. Up to t = TQL, the width of P is growing, since 
we proved (Ap^(i)) grows linearly over this time inter- 
val. Later on this width cannot decrease because of the 
locality of chaotic motion 15|. We assume t ^ Tgprcad- 
Then the width w of P is narrow enough for the spread 
of Sq to be negligible over a time t ~ w/km- Therefore 
{exp[ikjnSq{T\po,q{t"),t")])^ ~ 1 in the part of the in- 
tegration domain over t where P takes appreciable val- 
ues in (|l^, and B = limt^oo Z]m=i /o ^(^™ ~ 
Po,t"bo)dt" = J^J2f^^^M^P(v^ - poX\Po)dt", 
where the inverse Fourier transform was provided by the 
integral over t. 

Now, if t is large enough for P to be almost constant 
over the range [vm-LjVm+L] for all m's, we approxi- 
mate Y^.^_j^D„i+j\AVm+j\/\Vm+L+l - Vm-L\ ^ £'qL 
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and substitute the sum over Vm by an integral : B = 
2 lo I DQLPip - Po, t"\po)dpdt" = 2DQi^t. 

The general term of (|l^) can be estimated by a sim- 
ilar calculation. A sequence of two Fourier transforms 
is again recovered. After the first one, averages of the 
kind {expi[kmSq{T\po,q{t"),t") + ip,n + efn])* are found. 
They vanish as the constraint p{t") = po + Sp leaves al- 
most free the average on any two phases, and since Sq 
is negligible for t small. Therefore only B contributes 
to (Ap2(i)) which thus grows in a quasilinear way. This 
ends our proof of the quasilinear estimate for asymptotic 
times. 

Note that the conditional probability P permits to use 
the knowledge of initial quasilinear diffusion for proving 
it over asymptotic times only because we proved before 
that TQL ^ Tsproad- In contrast with the initial non- 
chaotic quasilinear regime, the number of modes acting 
on the particle increases with t. This agrees with the fact 
that the orbit visits an increasing number of resonances 
when time increases. 

VI. CONCLUSION 

Thus we prove the quasilinear character of the diffusion 
for the motion of a particle in a spectrum of large ampli- 
tude longitudinal waves. Our technique can be adapted 
to systems with a slow dependence of the quasilinear dif- 
fusion coefficient on p. As many Hamiltonian systems 
may be locally reduced to case (0) Hi, this further ex- 
tends its range of applicability and shows that the uni- 
versality class of quasilinear diffusion is broad. It also 
provides insight for the case where particles and waves 
are self-consistently coupled pTf . 

Higher order moments of Ap could be computed using 
a similar technique. Indeed, preliminary calculations in- 
dicate that the use of conditional probabilities should en- 
able one to retain after Fourier transforms the same terms 
for the moment of order k as in the case where q{t) is 
weakly dependent on any phase provided that k <C B~^, 
which yields a gaussian estimate. Proving the Gaussian- 
ity of / would also lead to a Fokker-Planck-Smoluchowski 
evolution equation for /. 

The value of B (which depends only on local aspects 
of the spectrum : A, fc, Sv) determines the time scale 
over which the quasilinear approximation holds. Given 
B 1, this time scale is t ^ tql. On the other 
hand, we require that the motion remains away from the 
boundaries Pmin and Pmax of the wave spectrum. Given 
the scaling (Ap^) ~ 2Dt, the boundary is reached for 
tbound ^ D-^APAv"^ - NPBtqi^. As M is independent 
of B^ one may let M ^ cx3 to ensure tbound to be as large 
as desirable. 

Comments by D. Benisti and A. Henriet on this work 
are gratefully acknowledged. YE thanks the organizers 
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